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Abstract 

We describe how to calculate the parton distributions f a (x,k 2 ,/j, 2 ), unintegrated over 
the parton transverse momentum kt, from auxiliary functions h a (x,kf), which satisfy 
single-scale evolution equations. The formalism embodies both DGLAP and BFKL con- 
tributions, and accounts for the angular ordering which comes from coherence effects in 
gluon emission. We check that the unintegrated distributions give the measured values of 
the deep inelastic structure function F2(x,Q 2 ). 

1 Introduction 



Conventionally deep inelastic lepton-proton scattering is described in terms of scale-dependent 
parton distributions, a(x, /i 2 ), where a = xg or xq. These distributions correspond to the density 
of partons in the proton with longitudinal momentum fraction x, integrated over transverse 
momentum up to k t = fi. They satisfy DGLAP evolution in /i 2 . The kinematic region k t < fi 
gives the leading In fj, 2 approximation to deep inelastic scattering. 

For less inclusive processes it is, however, necessary to consider distributions unintegrated 
over the transverse momentum k t of the parton. The unintegrated distributions have the 
advantage that they exactly correspond to the quantity which enters the Feynman diagrams 
and therefore allow for the true kinematics of the process even at leading order (LO). These 
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distributions f a (x,k^, /i 2 ) depend on two hard scalesQ k t and the scale \i of the probe. The 
scale /i plays a dual role. On the one hand it acts as the factorization scale, while on the other 
hand it controls the angular ordering of the partons emitted in the evolution [|TJ. 

Clearly it is desirable to also include ln(l/x) BFKL-type contributions in the evolution. Re- 
call that both DGLAP and BFKL evolution are essentially equivalent to ordered evolution in 
the angles of the emitted partons[|. In the DGLAP collinear approximation the angle increases 
due to the growth of k t , while in BFKL the angle {9 ~ k t /ke) grows due to the decreasing 
longitudinal momentum fraction as we proceed along the emission chain from the proton. The 
factorization scale \i separates the partons associated with emission from different parts of the 
process, that is from the beam and target protons (in pp collisions) and from the hard subpro- 
cess. For example it separates emissions from the beam (with polar angle 9 < 90°) from those 
from the target (with 9 > 90°), and from the intermediate partons from the hard subprocess. 
This separation was proved in Ref. [0] and originates from the destructive interference of the 
different emission amplitudes in the angular boundary regions. If the longitudinal momentum 
fraction is fixed by the hard subprocess, then the limits on the angles can be expressed in terms 
of a factorization scale fi which corresponds to the upper limitQ on the allowed values of the 
(s-channel) parton k t . 

Since the parton distributions depend on two scales we potentially have to deal with com- 
plicated (CCFM HI) evolution equations for the f a (x, fc 2 ,// 2 ) functions. Of course it is possible 
to work with two-scale distributions, but this is much more complicated and up to now 
has only proved practical with Monte Carlo generators ||. However, the evolution process is 
essentially controlled by one quantity, the emission angle, and on this basis we may expect to 
be able to obtain the distributions f a (x, kf, /i 2 ) from single-scale evolution equations. Therefore 
it should be possible to follow an analytic approach where the physical assumptions are much 
more evident and where, in principle, NLO corrections can be included. Moreover, in practice, 
it is much easier to use the same unintegrated distributions to describe different hard processes 
and to perform global parton analyses. 

The outline of this paper is as follows. The key observation is that the /i dependence of the 
unintegrated distributions enters at the last step of the evolution, and so we may use single- 
scale evolution equations. The procedure is first described in Section 2 in the case of pure 
DGLAP evolution, and then extended to include ln(l/x) effects in Section 3. In the latter case 
we use the solution of a single-scale equation which unifies DGLAP and BFKL evolution ||, 
and perform a final evolution step which brings in the dependence on the second scale. Ref. || 
also generated the two-scale unintegrated gluon from the same unified evolution equation, but 

This property is hidden in the conventional distributions as k t is integrated up to the scale /x. 

2 At LO we have strong ordering of the emission angles, . . .6^ <C 9i + i . . .; on the other hand if, at one step of 
the evolution 8i ~ then this contribution is included inside the NLO splitting function. 

3 The i-channel parton may have kt up to (i/z, characteristic of BFKL effects, whereas for LO DGLAP the 
s and i-channel partons are both limited by k t < \l. Of course, some k t > /i contribution will arise from the 
NLO splitting functions. 
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with a different procedure^ The unintegrated gluons obtained using the procedures described 
in Sections 2, 3 and Ref. || are compared in Section 4. In Section 5 we describe how the 
structure function F 2 is calculated from the unintegrated parton distributions, and in Section 6 
we discuss the relationship between the unintegrated and integrated distributions. Finally in 
Section 7 we give our conclusions. 



2 Unintegrated DGLAP partons 



It is informative to review how unintegrated distributions f a {x, k 2 ,fi 2 ) may be calculated from 
the conventional (integrated) parton densities, a(x,fi 2 ), in the case of pure DGLAP evolution. 
The procedure was explained in Ref. ||. We start from the DGLAP equation^ 

da(x, /i 2 ) 



d In /i 2 
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2tt 
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where in the first term a sum over all possible parent partons a 1 is implied. This first term on the 
right-hand-side describes the number of partons Sa emitted in the interval fj 2 < k 2 < fi 2 + 5fi 2 . 
Such emission clearly changes the transverse momentum k t of the evolving parton. If we were 
to neglect the virtual contribution in ([]]), then the unintegrated parton density would be given 
simply by 



da(x, /i 2 



<91n /i 2 

2tt 



J Paa'(z) a (z>kfj dz. 



(2) 



The virtual contribution in ([I]) does not change the parton k t and may be resummed to give 
the survival probability T a that parton a with transverse momentum k t remains untouched in 
the evolution up to the factorization scale. The survival probability is given by 



/, x ( a s(K 2 ) dk' 2 _ /-i-A A 



(3) 



a la Sudakov form factor. Thus the probability to find parton a with transverse momentum k t 
(which initiates a hard subprocess with factorization scale fi) is 



f a (x,kt,fi 2 ) = T a (kt,n) 



a s {k 



2tt 
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[z) a' 



- h 2 



dz 



(4) 



4 In this work we impose the angular ordering constraint in both the BFKL and DGLAP terms, and as 
a result do not have an exact equality between the integral up to /i 2 of the unintegrated distributions and 
the value of the integrated distribution. Ref. [^| takes the opposite approach; that is, exact equality with the 
integrated distribution is imposed and as a result angular ordering of the BFKL contribution is not complete. 
The difference is a NLO effect. 

5 For the g — + gg splitting we have to insert a factor z' in front of P gg (z') in the last integral of (Q) to account 
for the identity of the produced gluons. 
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It is at this last step of the evolution that the unintegrated distribution becomes dependent on 
the two scales, k 2 and fi 2 . 

We now have to take care to specify the value of the infrared cut-off A, which is introduced 
to protect the 1/(1 — z) singularity in the splitting functions arising from soft gluon emission. In 
the original DGLAP equation, (|l|), which describes the evolution of the integrated distributions, 
this singularity is cancelled between the real emission and virtual contributions. However after 
the resummation of the virtual terms, the real soft gluon emission must be accounted for 
explicitly since it changes the k t of the parton. Thus we have to find the physically appropriate 
choice of the cut-off A to provide the angular ordering of the gluon emissions^. 

In Ref. the cut-off was taken to be A = k t j ft. As a consequence the two-scale unintegrated 
distributions f a (x,k 2 ,fi 2 ) of || vanish for k t > fi, in accordance with the DGLAP strong 
ordering in kt- However we can do better and impose the more correct angular ordering in the 
last step of the evolution. It was shown in Refs. |1], || that this leads to a constraint on the 
scale /i, namely 

9(0-0') \i > zk t /(l-z). (5) 
Thus the maximum allowed value of the integration variable z is 

Zmax i ; \y) 

and the corresponding cut-off A = k t /(fi + k t ). Of course the same A must be used both in 
the real emission integral in (^ and in the survival probability T in (§){] In fact we shall see 
that the imposition of angular ordering at the last step of the evolution leads to physically 
reasonable parton k t distributions which extend smoothly into the domain k t > ft. 



3 Inclusion of ln(l/cc) effects 

We wish to generalize the above method to include the leading ln(l/x) contributions. Clearly 
there can be different forms of the 'unified' evolution equation summing up the leading DGLAP 
and BFKL logarithms, where the ambiguity is at the subleading level. The aim is to find a 
good prescription which is not too complicated, but which can account for all the physically 
relevant kinematic effects just at LO level. In other words we seek an equation which sums up 
the major part of the subleading corrections in a LO framework. 

Let us consider, for the moment, just the gluon distribution. Recall that the unintegrated 
distribution f(x, k 2 , ft 2 ) depends on two scales. As in the pure DGLAP case of Section 2 we wish 
to work in terms of a single-scale evolution equation, and then to restore the scale ft, and the 
full kinematics, at the last step of the evolution. This is illustrated schematically in Figure [l|. 

6 Although the splitting functions P gq and P qg are not singular at z = 1 it is natural to use the same 
prescription for both the quark and the gluon distributions. 

7 In equation (||), A = k' t /(fi + k' t ) is the appropriate cut-off for z'. 
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Figure 1: An illustration of our procedure, in which the evolution of a single-scale unintegrated 
parton is followed by a final step of the ladder which introduces dependence on the second hard 
scale, ii. 



For an analysis which incorporates BFKL effects, the appropriate single-scale distribution is 
the auxiliary function 

2 d(xg(x,k 2 )) 



h(x,k t [ 



(7) 



d In k 2 

Note that h(x, k\ ) is precisely the function which satisfies the BFKL equation in the low x 
limit. 

Both BFKL and DGLAP evolution correspond to angular ordering of the emission angles, 
and are single-scale equations. At LO, strong angular ordering automatically comes either from 
strong ordering in x (z -C 1) for BFKL or from strong k t ordering (k' 2 -C kf) for DGLAP. For 
the unintegrated gluon, we face a problem when k t ~ \x in the DGLAP framework, and similarly 
we have a problem when z ~ 1 for BFKL. Following the procedure of Ref. ||, we first neglect 
the subleading k t ~ fi and z ~ 1 effects to obtain and solve a unified BFKL/DGLAP equation 
for h(x, k 2 ). Then, in the last step of the evolution, we take account of the precise kinematics 
so that the kt ~ // and z ~ 1 domains are treated correctly. 

The unified equation for h g , which closely follows that presented in Ref. [Q, takes the form 
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where k' t = k t + (1 — z)q, see Fig. 
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Figure 2: part of the evolution chain. We commonly write fct for /c tn and then the parent's 
transverse momentum as k' t . The radiated transverse momentum is qt- Unified evolution is 
naturally performed || in terms of the rescaled transverse momentum q n = q tn /(l — z n ). 

We have introduced single-scale unintegrated quark auxiliary functions h q (x,fi 2 ) on the 
same footing as h g , 

V*,//) = »f», (9) 



and in (||) we sum over all 2np active flavours q of quarks and antiquarks with m q < /i. The 
h q distributions satisfy the equation 

h q ( x ,^) = h° q (x,f) + ^1 £ dz jT" |f [e(z-x) [p qg {z)h g {^-x) 



P qq {z)h q [-kf 



P qq (z)h q (x,k?)\, (10) 



which is the usual DGLAP equation for quark evolution written in terms of h q of ([|). The last 
term of (§) is the BFKL contribution which sums up all the leading (aglnl/a;)™ terms, while 
the remaining terms on the right-hand-side describe the conventional DGLAP evolution of the 
gluon distribution with respect to scale /x 2 . To avoid double counting, we have excluded the 
singular part of the P gg splitting function in the real emission DGLAP term and used 

2 N 

P(z) = P gg (z) - — °. (11) 

The 2Nc/z term is already included in the BFKL contribution to (|8|). 
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The driving terms, h°, which describe the low k 2 < k$ domain are given by 



K(x^ 2 ) = £ dz[Q(z-c 

+ P gg( z ) k o) J - zP gg (z)xg(x i k ) - P gq (z) Y,xq(x,k[ 



(12) 



and similarly for h®. The integrated input distributions a(x, k$) are not known and, as usual, 
must be determined from the data or from some non-perturbative QCD model. 

We emphasize that in the unified BFKL/DGLAP equation we choose the scale \x to be k t 
for the DGLAP contribution, to be consistent with the BFKL term (which is independent of \i 
at LO). Recall that angular ordering leads to a redefinition of the scale |5|] 

6(0-0') => n > zk t /(l-z). 

If this modified scale were to be adopted then it would be impossible to obtain a simple one- 
scale unified evolution equation throughout the whole x, k 2 domain. On the other hand, within 
the LO framework, we may omit the z dependence in the scale, so that the DGLAP part of 
the evolution becomes ordered in transverse momenta. We stress again that it is sufficient to 
implement the precise constraints coming from angular ordering (relevant to the z ~ and 1 
domains for the BFKL and DGLAP terms respectively) at the last step of the evolution. 

An advantage of the single-scale unified BFKL/DGLAP equation is that it is straightforward 
to incorporate a major (kinematical) part of the subleading order ln(l/a;) (BFKL) effects^ 
by imposing a consistency condition to ensure that the virtuality of the exchanged gluon is 
dominated by its transverse momentum squared f?J. This is achieved by the inclusion of the 
theta function 0(/i 2 — zq 2 ) in the real emission contribution shown in the last term of Note 
that other subleading effects arising from using the complete DGLAP splitting function and 
running as are automatically included in the unified equation. 

We see that the evolution equations for the auxiliary distributions h a (x, k 2 ) depend on the 
single scale kf. The dependence of the unintegrated distributions f a (x,k 2 , jJ 2 ) on the second 
scale n will enter when we consider the last step of the evolution. It is sufficient to ensure that 
the final emitted parton explicitly satisfies the requirements of angular ordering. The angular 
ordering conditions of the previous steps of the evolution are automatically ensured at LO by 
virtue of either the strong ordering in kt (in the DGLAP part) or the strong ordering in z (in 
the BFKL part). 

To ensure angular ordering in the last step of the evolution, we note that z is limited by 
d^). This condition implies z < /z/(/i + k t ), and so we take this as the upper limit of the z 

8 The large NLO BFKL corrections, which have recently been computed ||, appeared to have put the 
application of the BFKL framework into question. However a major part of the corrections is kinematic in 
origin and, when summed to all orders H using the theta function 0(/i 2 — zq 2 ), brings the BFKL approach 
back under control, see also @, [l(J . 
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integration in flTBD and (|T3D below. Thus the number of gluons produced at the last step (with 
transverse momentum k t which initiate a hard subprocess with factorization scale /i) is£] 



f g (x,k 2 ,fi 2 ) = T g (kt,n) 




- h' 2 
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(13) 



and the number of a particular quark species q, produced at the last step, is 



f q (x, k t , fx 
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This last step of the evolution is shown schematically in Fig. [T]. The unintegrated distributions 
f a {x, k 2 , fi 2 ) represent the number of partons 6a emitted in the S In k 2 interval from k\ to k 2 +5k 2 , 
and include the factors T a (k t , fj) of @, which give the probabilities that partons a = g,q with 
transverse momentum k t remain untouched in the DGLAP evolution up to the factorization 
(probe) scale \i. These survival probabilities T g and T q resum the virtual DGLAP contributions 
occurring in (§) and (|T0| ) respectively. 

It is important to note that the function h g (x,k 2 ) already includes the leading ln(l/x) 
virtual corrections, which have the effect of reggeizing the exchanged gluon; thus k t is the total 
momentum transferred via the Regge gluon trajectory. Moreover, the distributions f a (x, k 2 , /x 2 ), 
evolved in the final step from the auxiliary functions h a (x,k 2 ), also incorporate the virtual 



DGLAP contributions, via the survival probabilities T a (kt,n). The final expressions, ( |T3| ) and 
([H]), are thus more symmetric in that all the LO virtual corrections are included. That is, 
from a Feynman diagram viewpoint, the function f a (x,k 2 , /i 2 ) corresponds to the propagator 
of a t-channel parton of transverse momentum k t , initiating a hard subprocess at scale fi, in 
which all the LO virtual corrections to the parton distribution have been taken into account. 



4 The unintegrated gluon distribution 

We have described how to obtain the (two-scale) unintegrated parton distributions f(x, k 2 ,fi 2 ) 
from the solution h(x, kf) of a one-scale equation which unifies DGLAP and BFKL evolution. 
The link is (0), (|14]), which represent the last step of the evolution. Only at this stage does 
the scale /i of the subprocess, initiated by f(x, k 2 ,fi 2 ), enter. 

9 The low kt < ko domain of the integrals in ([l3]) , Jl4| ) should be treated as the driving terms h° in (|lj) . 
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It is informative to compare the unintegrated gluon distribution f g (x,k 2 , /z 2 ) with the be- 
haviour of the auxiliary function h g (x,k 2 ). Consider first the pure BFKL limit, in which 
DGLAP evolution is neglected, that is 0:5 In /i 2 <C 1, but a^lnl/x > 1. Then, in flTB|), the 
survival probability T g (kt,fi) = 1 and there are no DGLAP contributions. Thus we obtain 

f g (x,k 2 ,^ 2 ) = h g (x,k 2 ) = d^fop . (15) 



This is an expression which is frequently used at low x. 
When x is sufficiently large, the derivative 



d In k 2 1 DJ 

becomes negative, even for k t < \x. The reason is that the negative virtual DGLAP term exceeds 
the real emission DGLAP contribution, which is suppressed by the large lower limit z > x in 



( p~3| ) and (14). After the virtual contributions are resummed into the T a factors the 
unintegrated parton distributions remain positive everywhere. 

We already stated in Section 2 that imposing angular ordering on the last step of "DGLAP" 
evolution caused the distribution to extend smoothly into the k t > \i domain. The gluons 
also populate this domain due to BFKL diffusion in ln/c 2 . This raises the question of how 
much enhancement we find in this domain due to the inclusion of the BFKL contributions. 
To investigate this point, we compute the unintegrated gluon distribution f g (x,k 2 ,fi 2 ) in two 
different ways. 



(a) In the first approach, the function h g (x, k 2 ) of Ref. M is used as the auxiliary func- 
tion to drive the last-step evolution in ([□]) and flT4]). This incorporates essentially 
maximal BFKL effects.^ 

(b) In the second case, we use unintegrated "DGLAP" partons to drive the two-scale 
unintegrated f a via (f|). This approach is essentially pure DGLAP, but with the 
crucial modification that the cut-off in (|3|) and (|]) is motivated by angular ordering. 
Here we use the MRST99 1TTJ set of partons as input. 



Ideally, we should refit to the deep inelastic and related scattering data and perform global 
unintegrated parton analyses^ in terms of f a (x, k 2 ,fi 2 ), since now we have added an extra last 
step to the evolution. However, such global analyses are beyond the scope of the present paper. 

The unintegrated gluon distribution f g (x, k 2 , /i 2 ), obtained by the two alternative procedures 
(a) and (b), is shown in Fig. |3] for /i = 10 GeV. We can compare the continuous curves 

10 In comparison with , the evolution of h q in Ref. Q contains some resummation of the BFKL-like leading 
(as In l/x) n terms. 

11 In Section 5, we describe the theoretical calculation of jpg, the most important observable for constraining 
parton sets, from the two-scale unintegrated f g and f q . 
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of approach (a) directly with the previous unintegrated gluon distribution (dashed curves) 
calculated in Ref. ||, which also used the same auxiliary function h g (x, of ]|] as input. For 
kf < [L the new results tend to lie above the previous determination ||, while for k t > fi they 
are increasingly smaller as k t increases. In the present work, the decrease at large k t arises 
from the restriction z < z max = fi/(fi + k t ), whereas the earlier calculation, based on eq. (23) 
of Ref. ||, omitted the limit z < z max on the BFKL contribution, which causes f g to increase 
in the large k t domain at small values of x. 

Finally we compare the DGLAP-like unintegrated gluon, obtained in approach (b), with 
that of approach (a) which embodied BFKL evolution. That is we compare the dotted with 
the continuous curves of Fig. [3[ It is interesting to note that, with the more precise angular- 
ordered (CCFM |1|) cut-off from (|), A = k t /(fi + k t ), the unintegrated gluon distribution (b) 
generated from (|j) using the pure DGLAP MRST99 [TT| partons turns out to be very similar 
to the result (a) obtained using the auxiliary function h g . This coincidence can be explained by 
the facts (i) that both the DGLAP partons and the analysis of Ref. Q (which yielded h g ) fit the 
deep inelastic data well, and (ii) that both these input are used with the same angular-ordered 
constraint, @. We draw the conclusion that the role of angular ordering in the last step of 
evolution is particularly important, even more so than BFKL effects in the HERA domain. 



5 2*2 calculated from the unintegrated partons 



To check the reliability of our unintegrated parton distributions, and to demonstrate how to 
use these distributions in calculations of observables, we compute the deep inelastic structure 
function F 2 . We wish to treat the unintegrated gluons and unintegrated quarks on an equal 
footing as input to the subprocess cross sections, so we explicitly separate gluon and (direct) 
quark contributions to F 2 . 

The gluon contributes to F 2 via the quark box and crossed-box diagrams of Fig. These 
generate, via the g — > qq splitting, a sea quark contribution S q to F 2 of the form |T2], [| 

12^ 



(17) 



with 



S q (x,Q 2 ) = g J ^ jf 1 dp J d\ t a s (^)f g ^kl^ e(i-~ 



p 2 + (1 - pf 

The denominator factors are 
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(19) 
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Figure 3: plots of the fcf-dependence of the unintegrated gluon f g (x,kf, /z 2 ) for various values 
of x, at \i — 10 GeV. The solid curves are our version (a) of f g from for comparison we 
show with dashed lines the unintegrated gluon from || (as in ||, the dashed lines have been 
smoothed in the transition region k t ~ fi). Also we plot our "DGLAP" unintegrated gluon (b) 
from (Q) in dotted lines, which with the correct angular ordering cut-off is very close to the new 
f g , especially at high x. 
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Figure 4: The quark box, and crossed-box, diagrams which mediate the contribution of the 
unintegrated gluon distribution f g (x/z,k 2 ,fi 2 ) to F 2 . 

We may exploit the symmetry of the integrand in fllTf ) under K t —>■ n t — k t and /? — > 1 — /3 to 
rewrite {. . .} as 

^ {[*■+<!-«'] (| - + H + W-fll - 5^)}- 

(20) 

The summation in ([i~7|) is over massless u, d, s quarks and a c quark of mass m c = 1.4 GeV; 
there is no need to sum w, d, s, c in addition, because as (|1|) is written, S u , say, is the 
contribution of a gluon via a u quark box of any momentum. The variable (3 is the light-cone 
fraction of the photon momentum carried by the internal quark. The variable z is the ratio of 
Bjorken x and the fraction of the proton momentum carried by the gluon. It is specified by the 
relation 

1 fa - (1 - (3)k t ? + m 2 q kl 

z + P(l-P)Q 2 Q 2 

(21) 

+ m 2 kj + k 2 - 2Kt -kt + ml 
+ (1-(3)Q 2 (3Q 2 

which is obtained by requiring the outgoing quarks to be on-shell. Following Ref. M, we choose 
the scale \i which controls the unintegrated gluon distribution and the QCD coupling as to be 

/i 2 = k\ + k 2 + m 2 . (22) 

Care is needed in separating this calculation into perturbative and non-perturbative regions. 
We impose a cut-off k t > k for a legitimate perturbative calculation of (Pf) . The smallest cut- 
off we can choose is the minimum (that is, initial) scale of the function h g (x,k 2 ) from which 
our two-scale distributions derive. Thus ko is of order 1 GeV. For the contribution from the 
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region of k t < k we approximate 



■o dk 



TT J< 



g\X, k t , jj, ) 



remainder 



xg(x, kl) T g (k ,/j) 



(23) 

The dependence on 



where a can be taken to be any value representative of the interval (0, k 
the choice of a is numerically unimportant. 

Now we have to add the direct quark contributions to F 2 , which come from the unintegrated 
quark distributions f q (x, k 2 , fi 2 ). If a quark, initially with x/z and perturbative transverse 
momentum k' t > k , splits to a radiated gluon and a ' final Q quark with Bjorken x and transverse 
momentum K t , then this final quark can couple to the photon and contribute to F 2 as: 



,-i?(pert) 



(x,Q 2 
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dn 2 a s (n 2 ) 



2tt 



dk 2 



QKQ+kt) 
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fg ( ~-> k't > Qj + fg k 2 , Q 2 



(24) 



where here we have written the antiquark contribution explicitly. As in (P^), the upper limit 
of the z integration reflects the angular-ordered constraint of ([|) during the quark evolution. 

Again we need to account for the non-perturbative domain k t < ko- The initial (integrated) 
quark distribution xq(x, k 2 ) drives our final contribution. Physically the only remaining dia- 
grams that we have not included are those in which a quark (or antiquark) from this initial dis- 
tribution does not experience real splitting in the perturbative domain, but interacts unchanged 
with the photon at scale Q. Hence we write a Sudakov-like factor T q (ko,Q) to represent the 
probability of evolution from ko to Q without radiation. 



-,<jr(non— pert) 



x,Q 2 ) 



(xq(x,kl) + xq(x,klfj T q (k Q ,Q). 



(25) 



To avoid double counting, it is important to put a lower limit on K t in both (^) and (p4]) 



by enforcing Q(k 2 — k$) in the integrations. Without this lower limit on the final transverse 



momentum K tl equations (|18D and (^J) would partially include low transverse momentum K t 
quark contributions which are best incorporated in (pB"D, whether they originate from partons 
with k t > ko or not. 

The structure function F 2 (x, Q 2 ) is given by the sum of the gluon-initiated contribution (|17|). 
and the quark terms, ( f24| ) and (p5|). In Fig. || sample results^ shown by the continuous curves, 
are compared with deep-inelastic structure function data. The gluon and quark components are 
shown by the dashed and dotted curves respectively. As expected the dominant contribution 
at small x comes from the unintegrated gluon via the quark box and crossed-box contributions, 
whereas at large x the quark terms dominate. 

12 By 'final' we mean the struck quark just after the P qq splitting. 

13 We show in Fig. || results obtained using the unintegrated distributions evaluated from (^) with the MRST99 
partons @, that is, version (b) of f g and f q as discussed in Section 4. 
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We emphasize that, in the present work, the curves for F 2 are not the result of a fit to the 
structure function data. Rather they have been obtained by using single-scale functions, origi- 
nally fitted more directly to F 2 data, as plausible input to our 'last-step' evolution procedure, 
which generates two-scale unintegrated distributions f a (x, k 2 , /i 2 ). Then we use the uninte- 
grated distributions to compute F 2 via (|17|), (^4j) and For the inclusive observable F 2 , we 
would expect that the insertion of this extra evolution step would not appreciably disturb the 
description, since it essentially redistributes the distributions in k t space. We see from Fig. |] 
that indeed this is the case. 



6 Relation of / to integrated partons 

It is important to scrutinise the relationship between the new unintegrated partons f a (x, k 2 , /i 2 ) 
and the conventional integrated parton distributions a(x,/i 2 ), as obtained in global analyses 



such as |pj{| . First we emphasize that we may use either the integrated distributions or the 
unintegrated distributions to describe both inclusive (like F 2 ) and exclusive processes. The 
framework based on the unintegrated distributions is a bit more complicated. However it 
accounts for the precise kinematics of the process and an important part of the virtual loop 
corrections, via the survival factor T, even at LO. On the other hand, if we work with integrated 
partons we have to include NLO (and sometimes NNLO) contributions to account for these ef- 
fects. These differences appear to cause a discrepancy between the integrated and unintegrated 
approaches. As we explain below, this is to be expected since it arises from simplifications of 
the LO formalism due to the neglect of terms which are moved into the NLO contribution. 

An important equation, sometimes cited as the defining property of unintegrated partons 



is 

, ^ 2 dk 2 

a(x, jj, 



f ^ fa f), (26) 



where a represents xg or xq. This is in fact the first equation of Ref. ||. In the BFKL limit, 
the \i dependence of / vanishes and we have f g (x, k 2 ,fi 2 ) — > h g (x, k 2 ) as in fll5l). In this case, 
( |26| ) is clearly satisfied. However, in general the situation is complicated by the two separate 
momentum scales k t and fi. The unintegrated partons f g of Ref. || were explicitly constructed 
to have the property fl2"6|) , in the sense that the integral of f g over the transverse momentum up 
to the scale fi would be the same as the integral of the input auxiliary function h g (x, k 2 ) up to 
the same scale. In contrast, numerical integration over k t of the new unintegrated partons f g 
and f g presented in this paper (both versions (a) and (b) of Section 4) shows that ([ZED i s om y 



approximately true.0 We typically find a discrepancy of order 25% between the right-hand 
side of fl2"6|) and the single-scale distribution that has been used to generate f a . 



14 Note that we cannot compute ( Pq ) as it is written, because we can only define the unintegrated function in 
the regime of perturbative k t > kg. The comparison that is made is between the integral from k^ to /i 2 and the 
quantity a(x, /i 2 ) — a(x, fcg). 
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Figure 5: This is not a fit but the results of using our "DGLAP" unintegrated partons (b) 
to calculate F 2 ; the gluon-originated contributions are shown as dashed lines and the quark- 
originated parts are shown as dotted lines. Recent data are plotted fl^[| , and compare well with 
the sum of the gluon and quark contributions (solid curves), especially at high Q 2 . 
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In order to eliminate the discrepancy we may adjust the upper limit /i 2 of the integral 
in fl26|) to c 2 /i 2 . The introduction of c is equivalent to a NLO correction for the integrated 
partons. Typically in the "DGLAP" case (approach (b) of Section 4) we require c = 0.6 — 0.8 
to reproduce^ the original MRST integrated gluon in the domain \i — 5 — 10 GeV and x < 0.01 
(and c ~ 0.4 for approach (a), which embodies BFKL effects). A value c < 1 compensates for 
the over- large virtual loop correction included in the integrated DGLAP partons on account of 
the absence of the cut-off A in the integral over z'. At smaller values of x the corrections due 
to the cut-off are mainly cancelled between the virtual and real DGLAP contributions. As x 
increases, the virtual contribution (the second term on the right-hand-side of ([!])) increasingly 
dominates and we have to choose smaller values of c. Eventually in the domain x > 0.1, 
~ 10 GeV the main contribution comes from the input and to compensate the absence of the 
z cut-off in the conventional form of the DGLAP equation we must change the input itself. 

To summarize, the discrepancy between the integral (p6|) of the unintegrated parton function 
and the original integrated distribution is not a cause for concern. Conceptually, there are two 
different roles for single-scale distributions in the description of data for inclusive observables 
(where partonic transverse momentum is integrated out). The first role is the traditional one, in 
the framework of collinear factorization, whereby integrated parton distribution functions are 
fitted directly to the data. The second role is demonstrated in this paper (following || and ||), 
where we use single-scale functions as input to the last-step procedure, see for example (|J). We 
have been forced to introduce a new formalism for calculating F 2 (x, Q 2 ), in which unintegrated 
parton functions are understood to be the fundamental objects; we have emphasized the need 
to perform a new global fit to data in terms of the new functions f a . After this, we do not 
expect the input single-scale function a on the left-hand side of (|26|) to equal the integral of 
f a up to /i 2 , since a itself is not fitted directly to the data, but rather is used as input for the 
last step of the evolution, which embodies a crucial angular-ordering constraint unique to this 
last step. Thus our single-scale or 'auxiliary' function is not a traditional parton distribution 
function, but simply an intermediate function. 



7 Conclusions 

Parton distributions, f a (x,k 2 , fi 2 ) unintegrated over the parton k t are the basic quantities for 
describing processes initiated by hadrons. An essential ingredient in this description is the 
existence of the kt factorization theorem WM. The unintegrated distributions depend on two 



hard scales — the transverse momentum k t and the factorization scale \i. The scale /i drives the 
angular ordering during the evolution which arises from the coherence of the gluon emissions. 

Here we develop a new formalism to determine the unintegrated parton distributions, 
f a (x, k 2 , /i 2 ), which embodies both the leading \xiQ 2 (DGLAP) and lnl/x (BFKL) effects, 
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This identification is only approximate since the last step of the evolution is not included in the comparison. 
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as well as including a major part of the sub-leading Inl/x contributions. An important ob- 
servation is that, at leading order, the two-scale functions f a (x, k 2 ,fi 2 ) may be calculated from 
auxiliary functions h a (x,k 2 ) which satisfy single-scale evolution equations, since the angular 
ordering restrictions, controlled by p,, become important only at the last step of the evolution. 
The equation for the auxiliary function h g (x,k 2 ) was formulated, and the distributions fitted 
to the data, in Ref . || . These single-scale equations were also devised to include all the leading 
ashiQ 2 and a^lnl/x contributions, and a major part of the sub-leading hxl/x effects. 

In other words, the 'unified' evolution equations for h must be supplemented by a final 
evolution step in which the p dependence of the unintegrated parton distributions enters via 
the angular-ordering constraint. The situation is summarised diagrammatically in Fig. [I]. The 
procedure offers a considerable simplification in the determination of physically realistic un- 
integrated parton distributions f a (x,k 2 , p 2 ) over the full x, p 2 perturbative domain, including 
the true kinematics even at leading order. As expected, the gluon and sea quark distributions 
extend into the k t > p region more and more as x decreases. We have compared the new 
unintegrated distributions with those given by previous prescriptions [|, ||. As compared to 
Ref. ||, the new formalism gives a consistent treatment of angular ordering, which leads to 
the imposition of the z integration limit z < p/(p + k t ). As a consequence the distributions 
f a (x, k 2 , n 2 ) decrease faster than those of || for large k t , particularly at small x. An interesting 
result is that the unintegrated distributions obtained via h a (x, k 2 ) of M are not very different 
from those obtained via (|J) using conventional DGLAP partons — compare the continuous and 
dotted curves in Fig. |. It thus appears that the imposition of the angular-ordering constraint 
is more important than the BFKL effects. This observation has the practical consequence that 
reasonably accurate predictions for observables can be made using the much simpler, though 
less complete, prescription of @. 

Finally, we used the new unintegrated distributions to calculate the deep-inelastic structure 
function F%. We also showed the gluon- initiated and quark contributions separately, which, as 
expected, dominate at small x and large x respectively. Recall, from J|, that the rise of the 
gluon at small x is driven by perturbative QCD, which was assumed to have a non-perturbative 
input which is 'flat' in x. We emphasize that we have not fitted to the deep inelastic data. 
Nevertheless, Fig. [5] shows that the existing distributions give an adequate description, and 
therefore they may be used to evaluate other hard processes, such as bb and large qt prompt 
photon production in high energy pp (or pp) collisions. It is important to use unintegrated 
distributions for such exclusive reactions. 
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